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XT) • ABSTRACT 

O^l i We quantify the importance of mass accretion during AGN phases in the growth of 

supermassive black holes (BH) by comparing the mass function of black holes in the 
local universe with that expected from AGN relics, which are black holes grown entirely 
with mass accretion during AGN phases. The local BH mass function (BHMF) is 
estimated by applying the well-known correlations between BH mass, bulge luminosity 
and stellar velocity dispersion to galaxy luminosity and velocity functions. We find that 
different correlations provide the same BHMF only if they have the same intrinsic 
dispersion. The density of supermassive black holes in the local universe which we 
[ estimate is pbh = ^i-^tiA K.i ^0 Mpc"^. The rehc BHMF is derived from the 

■ continuity equation with the only assumption that AGN activity is due to accretion 

onto massive BH's and that merging is not important. We find that the relic BHMF 
at z = is generated mainly at 2; < 3 where the major part of BH's growth takes 
place. Moreover, the BH growth is anti-hierarchical in the sense that smaller BH's 
O \ {Mbh < W'^Mq) grow at lower redshifts {z < 1) with respect to more massive one's 

{z ^ 1 — 3). Unlike previous work, we find that the BHMF of AGN relics is perfectly 



> 



o 

^ ■ consistent with the local BHMF indicating the local black holes were mainly grown 

' during AGN activity. This agreement is obtained while satisfying, at the same time, 

the constraints imposed from the X-ray background. The comparison between the local 
and relic BHMF's also suggests that the merging process is not important in shaping 
rS I the relic BHMF, at least at low redshifts (z < 3), and allows us to estimate the average 

■ radiative efficiency (e), the ratio between emitted and Eddington luminosity (A) and 

the average lifetime of active BH's. Our analysis thus suggests the following scenario: 
local black holes grew during AGN phases in which accreting matter was converted 
into radiation with efficiencies e — 0.04 — 0.16 and emitted at a fraction A = 0.1 — 1.7 
of the Eddington luminosity. The average total lifetime of these active phases ranges 
from ~ 4.5 x 10^ yr for Mbh < IO^Mq to ~ 1.5 x 10® yr for Mbh > IO^Mq but can 
become as large as ~ 10® yr for the lowest acceptable e and A values. 

Key words: black hole physics - galaxies: active - galaxies: evolution - galaxies: 
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1 INTRODUCTION 



Since the discovery of quasars JSchmiddllOel it has been 
suggested that Active Galactic Nuclei (AGN) are pow- 
ered by mass accretion onto a supermassive Bl ack Hole 
(BH) with mass in the range 10'^-10^° MfT. iSalpetej 
11964 ': 'Zel'do vich fc Novikovlll96l lLvnden-Bel^ll969^ . This 
paradigm combined with the observed evolution of AGNs 
implies that a significant number of galaxies in the local 
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universe should host a supermassive BH fe.g. ISoltanlll982l: 
ICavaliere fc Padovanilll989HChokshi fc Turnejll992h . 



BH's arc now detected in ~ 40 
gas and stellar dynamical methods 



Supermassive 
:alaxies through . 
Kormendv fc Gebhardjl200ll: iMerritt fc Ferraresdl2001bl) 



Some galaxies are qu iescent (e.g. M32. [ van^^^Mare^^ill 
[W9I; NGC 3250, iBarth et alJ l200j^ and 



some are 
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mildly or strongly active (e.g. M87, 
19971: IMacchetto et al.l Il997l: Centa.urus A , 
20011 Cygnus A. iTadhunter et alJ 120031) . The mass of 
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the BH correlates with some properties of the host 
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J alaxy such as spheroid (bulge) luminosity and mass 
Kormendv fc Richston j|l995l: iMagorrian et alJll998l) . light 
concentration (Graham et al. 200 1) and with the cen- 
tral stellar velocity d ispersion ( Ferrarese fc Merritti I2OOOI : 
iGebhardt et al.ll2000h . The latter correlat i on wa s thought 
to be the tightest but iMarconi fc Hun3 l)2003f) have re- 
cently shown that all the correlations have similar in- 
trinsic dispersion when conside ring only galaxies with 
secure BH detections (see a lso iMcLure fc Dxinlopl I2OO2I : 
lErwin. Graham, fc Caonll2003l) . Overall, the dispersion is of 
the order of 0.3 in logMeH at a given value of Lbui, Afbui or 
cr* . The existence of any correlations of BH and host galaxy 
bulge properties has important implications for theories of 
galaxy formation in general and bulge formation in partic- 
ular. Indeed, several attempts at explaining the origin of 
these correlations and the difficulties/constraints that they 
pose to gala xy formation models can be found i n the lit- 
eratu r e (e.g.lSilk fc R,eeslll998l: ICattaneo. Haehnelt. fc R.eei 



199c ; iHaehnelt fc KauffmannI 120001: ICiotti fc van Albada 
200 J : ICavaliere fc Vittorinil l2002l : lAdams et alJ I200I and 



references therein). 

To date, all newly found BH's have masses (or upper 
limits) in agreement with those expected from the above 
correlations suggesting that all galaxies host a massive BH 
in their nuclei. By applying the correlations between Mbh 
and host galaxy properties it is then possible to estimate the 
mass function of local BH's or, more sim ply, their total mass 
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Eooi 


lYu fc Tremaind l2002l: 


2002l:lAller & Richstond 


20021). 



" ^^m^esd 



It is important to verify if local BH's are exclusively 
relics of AGN activity or if other mechanisms, such as merg- 
ing, play a role. Hereafter we will call 'AGN relics', or simply 
relics, those black holes which grew up from small seeds (1- 
10^ Mq) following mass accretion during AGN phases. For 
instance, an AGN relic of 10^ Mq is different from a BH 
of the same mass but which was formed by the merging of 
many smaller BH's. 

A simple comparison between local and relic BH's was 
perfor med by Salucci et al. ( 19991) and iFabian fc Iwasawal 
l|l999|) who determined pbh from the observed X-ray back- 
ground emission. A revised est imate was obtained by 
lElvis. Risaliti. fc Zamoranil ^2002^ who found pBH = (7.5 — 



16.8) X 10^ (e/0.1)"' X lO'^ Mq Mpc 
the estimate from local black hole s ( 
lO'^ M ffi Mpc~^ lYu fc Tremaind 



in disagreement with 
3BH = (2.5 ±0.4) hi. 



20021: lAller fc Richstoiid 



I2OO2D . lElvis. Risaliti. fc Zamoranil ll2002ll thus suggested 
that, in order to reconcile this discrepancy, massive BH's 
should have large accretion efficiencies (i.e. larger than the 
canonically adopted value of e = 0.1), hence they should be 
ra pidly rotating. A mo r e deta iled comparison was performed 
bv lMarconi fc Salvatil ll2002ll who found an agreement be- 
tween the black hole mass functions (hereaft er BHMF) of lo- 
cal an d reli c black holes. Re cently, however. I Yu fc Tremainj 
J2002h and lFerrares3 i2002l) found a disagreement at large 
masses (Mbh ^ 10* M©) where more AGN relics are ex- 
pected relative to local BH's. As previously stated, a possi- 
bility to reconcile this discrepancy is to assume high accre- 
tion efficiencies but, clearly, this issue is still much debated. 

The relation between AGN relics and local BH's 
is also being studied in the framework of coeval evo- 
lution of BH and host galaxy. Several physical mod- 



els have been proposed in which the fueling of the 
BH, hence the AGN activity, is triggered by merging 
events (in the context of the hierarchical struct ure for- 
mation paradigm, see for instance ,Kau ffmann fc Haehneld 
2OOOI: IVolonteri. Haardt. fc Madaull2003l: iMenci et al.ll2003l: 
Wvithe fc Loebll2003l :l Hatziminaoglou et al .1 120031: iHaehneM 
2003h or is simply directly relate d to the star forma- 
tion history of the host gala xy (e.g. |D i Mattco et al."200j; 
iGranato et al] b004: Haiman. Ciotti. fc Ostriker 2003) . The 
BH has then a feedback on the host g alaxy through the 
energy released in the AGN phase (e.g. ISilk fc Reedll99'3 : 
Blandford 1999; Begelman 2003). As a result of this dou- 
ble interaction (galaxy feeding the BH - AGN feedback on 
the galaxy), these models can in general reproduce both 
the observed BH-host g alaxy correlations and the AGN lu- 
minosity functions (e.g. Haehnelt, Nataraian, fc ReeaM9g; 
Monaco, Salucci, fc Dancsc 2000; Nulscn fc Fabian! I2OO0I 
and previous references). However, this big effort in mod- 
eling cannot uniquely answer the question if local BH's are 
relics of AGN's, since a wide range of models with many 
different underlying assumptions cannot be ruled out with 
the available observational constraints. 



The aim of this paper is to investigate the assumption 
that massive black holes in nearby galaxies are relics of AGN 
activity by comparing the local BHMF with that of AGN 
relics. We remark that in this paper we do not build a phys- 
ical model of the coevolution between central BH and host 
galaxy but we compare differential and integrated mass den- 
sities (local BH's) with differential and integrated energy 
densities (AGN's), with the only assumption that AGN ac- 
tivity is caused by mass accretion onto the central BH. 



We refine the analysis bv lMarconi fc Salvatil i2002l) and 

evaluate the discrepancies found by other authors between 
local and relic BHMF's. In Section |5| we estimate the local 
BHMF by applying the known correlations between A^bh 
and host galaxy properties to the galaxy luminosity and 
stellar velocity dispersion function. We also check the self- 
consistency of the results, and show that different Mbh- 
host-galaxy-properties relations provide the same BHMF 
within the uncertainties. In Section |3 we use the continu- 
ity equation to estimate the BHMF of AGN relics and in 
Section 01 we compare local and relic BHMF's, and find 
that local BH's are consistent with AGN relics. We then 
show (Sec. 01 that the energetic constraints inferred from the 
X-ray background (XRB) are also satisfied and that there 
is no discrepancy between pbh of local BH's and that ex- 
pected from the XRB. In Sections |S| |7| and|H|we discuss 
constraints on the accretion efficiency (e) and on the Ed- 
dington ratio (A = L/LEdd, where L is the AGN luminosity 
and Z/Edd is the Eddington luminosity of the active BH) , and 
we estimate the accretion history and the average lifetime of 
massive BH's. We summarize our results and we draw our 
conclusions in Section |^ 



In this paper we adopt the current "standard" cosmo- 
logical model, with Hq = lOO/i kms"'^ Mpc~^ and h = 0.7, 
flM = 0.3 and = 0.7. 
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Table 1. Adopted galaxy properties per morphological types. 





E 


SO 


Sab 


Scd 


Notes 


Morph. Fraction 


0.11±0.03 


0.21±0.05 


0.43±0.07 


0.19±0.07 


1 


m(Bulge)-m(Total) 





0.64 ±0.30 


1.46 ± 0.56 


2.86 ±0.59 


2 


m(Bulge)-m (Total) 





0.60 ±0.50 


1.78 ± 1.01 


2.78 ±1.21 


3 



1. From Ipiikugita. Hoganj _^_Pccblcs 1998 (sec text for details). 

2. B band. Estimate bv TMl^fc Riclistoric 1.2002) based on data f rom lSimien &: de Vaucouleur j|l986l . 

3. H band. Based on data from lHunt. Pierini. &: Giovanardill2004i 



2 THE MASS FUNCTION OF LOCAL BLACK 
HOLES 

The first step of the analysis presented here consists in the 
determination of the mass function of local BH's, i.e. black 
holes residing in nearby galaxies. The sample of galaxies 
with dynamically measured BH masses is small (~ 40) and 
not selected with well defined criteria. Thus it is useless for 
a direct determination of the local black hole mass func- 
tion. However, the BHMF can be derived by applying the 
existing relations between Mbh and host galaxy properties 
to galaxy luminosity or velocity functions. After we have 
described the adopted formalism, we will verify the consis- 
tency of the Mbh-o"* and A/BH-ibui relations in providing 
the same BHMF within the uncertainties. We will then es- 
timate the BHMF for early and all galaxy types, showing 
that different galaxy luminosity /velocity functions provide 
BHMF's which are in agreement within the uncertainties. 

2.1 Formalism 

We describe here the simple formalism which is commonly 
used to derive the BHMF from galaxy luminosity or velocity 
functions. 

We define (j>{x)dx as the number of galaxies per unit 
comoving volume with observable x (e.g. luminosity L, or 
stellar velocity dispersion a) in the range x,x -\- dx. The 
observable y (e.g. Mbh, the BH mass) is related to x through 
the log- linear relation logj/ = a + 6 log a; and A(y) is the 
intrinsic dispersion in logy at constant log a:. Assuming a 
log-normal distribution then 



P(logy| logx) 



27rA, 



exp 



1 / log y — a — b log x 
'2 I A^) 



(1) 



where P(logy| \ogx)d\ogy is the probability that y is in the 
range log y, log y + d log y for a given log x. Thus the number 
of galaxies with x,y in the ranges x,x -[- dx and y,y + dy is 



^y, x) dx dy = £0£EA}£Slldy ^(^^^ 



(2) 



ylnlO 

${y)dy, the number of galaxies with y in the range y,y + dy, 
is thus the convolution of and P: 

P(logj/j logx) 



$(j/) = 



■ (t>{^) dx 



,0 2/ In 10 

In the limit of zero-intrinsic dispersion: 

$(y) = 



10" 6 



(3) 



(4) 



for instance, x with L, the spheroid luminosity, the total 
mass density in BH's is simply 



Pbh 



M${M) dM ■ 



dM 



dL- 



In 10 



P(logA/| log L)0(L) (5) 



'0 Jo 
and in the zero intrinsic dispersion case: 

(pbh)o = 10" / lV(^) dL (6) 

where M — Mbh for simplicity. 

Inverting the order of integration in Eq.j^and integrat- 
ing on Mbh one finally gets 

i[A(AlRH) In 101^ 

pBH = e2^ ^ bh; j p^^^ 



(7) 



A non-null A(AfBH) increases pbh by a factor 
exp[i(A(MBH) In 10)^] with respect to the zero dis- 
persio n case. 



a fact already noted by IYu fc Tremaind 



2.2 Consistency of the A/bh-ct* and MBH-^bui 
Relations 



a uthors (e.g. ISalucci et alll 1999HMarconi fc Salvatil 
'Yu & Tremaind l2002l : iFerraresd l2002l : 



with log y = a + b log x. After substituting y with Mbh and, 



Sever al 
'2002; 

AUer fc Richstone. .,2002i) adopted the method just de- 
scribed to determine the mass function of local BH's. In 
the most recent works the Mbh-o"* correlation has been 
preferred to M BH-J/bui on the ground that it is tighter 
and, moreover, IYu fc Tremaind ll2002h found a factor 2 
discrepancy between the values of pbh determined by 
applying the A/bh-o"* and AfsH-ibui relations. Thus, we 
first investigate this inconsistency of the AIbh-o'* and 
AfflH-ibui relations by examining the BHMF's derived by 
applying the two relations to the velocity and luminosity 
functions obtained from the same sample of galaxies. 

We consi der the SPSS sample of 9000 early type 
galaxies from iBernardi et al.l l2003al for which the lumi- 
nosity and velocity functions were determine d indepen- 
dently (IBernardi et al.ll2003bl: ISheth et alJl2003^ ■ By 'inde- 
pendently' we mean that the velocity function was derived 
by directly measuring a* of all galaxies and not by applying 
the Faber- Jackson relation (hereafter FJ) to the galaxy lumi- 
nosity function (e.g. Gonz alez c t al . 2000 ; AUcr fc Richst on3 
[2002; Fcrrarcsc 2002). S heth et'alj toO^ compared their ve- 
locity function with that derived applying the FJ relation to 
the luminosity function and found that it is fundamental to 
take into account the intrinsic dispersion of the relation in 
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Early Type Galaxies Early Type Galaxies 




7.0 7.5 8.0 8.5 9.0 9.5 10.0 7.0 7.5 8.0 8.5 9.0 9.5 10.0 

Log Mbh [M©] Log Mbh [M©] 



Figure 1. (a) Local BHMF for early type galaxies based on the SDSS sample of lBernardi et aT] j20033)- The shaded area and error bars 
("barred") indicate the 16 and 84% percentiles of 1000 Montecarlo realizations of the local BHMF and correspond to ±lcr uncertainties 
(see text). The A's indicate the assumed intrinsic dispersions of the Mrh-Cj- or MRH-^hni relations, (b) Local BHMF for early type 
galaxies derived using the luminosity functions from different surveys iBernardi et al.ll2003tjl iMarzke et alj|l994l: iKochanek et alJl200ll: 
iNakamura et allbOO.'j) and the A^BH-^bul relation with A(Mbh)=0.3. The shaded area and errorbars have the same meaning as in (a). 



order to obtain the correct velocity function. It is conse- 
quently expected that the same might apply to the determi- 
nation of the BHMF. 

To derive the BHMF using the 71/bh-(t* relation we ap- 
ply Eq. 13 For the galaxy velocity function (j) we use the 
param eterized form for early type galaxies bv ISheth et alJ 
i2003h . For the Mrh-c* relation we use t he coefficients ob- 
taine d by fitting the 'Group 1' galaxies of lMarconi fc Hund 
l2003l (i.e. galaxies with 'secure' BH mass determinations): 

log Mbh = (8.30 ± 0.07) + (4.11 ± 0.33)(log a, - 2.3) (8) 

The slope is in agreement with that of iTremaine et alJ 
|2002| (4.02 ± 0.32) but there is a larger normalization 
(8.13 ± 0.06 i.e. ~ 0.2 in logA/Bn). This choice of the co- 
efficients for the A/bh-o"* relation is made for consistency 
since, in the following, we use the coefficients of the Mbh- 
Lbui relation determine d for the same 'Group 1' galaxies by 
iMarconi fc HuntH2003h . The local BHMF's derived with the 
AfBH-c* relation assuming intrinsic dispersions A(AfBH)=0 
and A(AfBH)=0.30 are shown in Fig. The shaded area 
and errorbars indicate the 16 and 84% percentiles of 1000 
Montecarlo realizations of the local BHMF. The realizations 
were obtained by randomly varying the input parameters 
assuming that they are normally distributed with la uncer- 
tainties given by their measurement errors. For the velocity 
function we have considered only the error on cf)* , the num- 
ber density at a* , since th e other errors are strongly cor- 
related among themselves JSheth et alj|2003l) . The 16 and 
84% percentiles indicate the ilcr uncertainties on the loga- 
rithm of the local BHMF, whose values from the Montecarlo 
realizations are normally distributed at a given Mbh- 

To derive the BHMF using the MBH-ibui relation we 
again apply Eq. |21 The galaxy luminosity function (j> by 



IBernardi et"al] ^2003b^ is given as a function of the total 
galaxy light. Since the correlation of Mbh is with bulge 
light, we need to apply a correction Am in the case of SO 
gal axies to transform from total to bulge luminosity. Follow- 
ing '^^T'^T'TVemAin^ (|20o3), the luminosity function for the 
bulges of SO galaxies is directly given by 

/"so 

(j)b{m) = J, ' , , <t>{'m - Am) (9) 
./e + ,/so 

where (j) is the luminosity function of early type galaxies, m 
is the absolute bulge magnitude. Am = m — mtotai and /e, 
/so are the fractions of E and SO galaxies with respect to 
the total galaxy population. The galaxy type fractions and 
Am used in this paper are sho wn in Table [J The morpho- 
logical type fractions are from iFukueita. Hoean. fc Feeble^ 
119981 (their Sbc fraction has been evenly split between Sab 
and Scd). The la uncertainties are a conservative estimate 
we made after comparing various determinations of the mor- 
phological ty pe fractions available in the literatu re (see the 
discussion in iFukueita. Hoean. fc Feeble j [iggSl . T he Am 
for th e B band are th ose estimated bv lAl ler fc R ichstond 
J2002h by rebinning the|S imien fc de Vaucouleu^^^S^ 'data 
in the appropriate bins of galaxy types. The Am for the H 
band are based on data from ,Hunt. Picrini, fc Giovan ar'dil 
I2OO4I and t he values for SO galaxies were taken from the 
analysis bv lMarconi fc Hunt! i2003l) . The Am values are lit- 
tle dependent on the photometrical band at least within 
the considerable scatter, therefore in the following analy- 
sis we will always use the B band data regardless of the 
photometric band in which (j) -was measured. Applying the 
above c orrections to the early-ty pe r* band luminosity func- 
tion by IBernardi et alj i2003lJ) and the color transforma- 
tion r* — K = 2.8±0.2 we obtain the luminosity func- 
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Table 2. Mass density of local BH's for early and all galaxy types (Col. 3 and 4, respec- 
tively). Col. 1 indicates the galaxy luminosity or velocity function which was combined with the 
A/eH^host-galaxy-property in Col. 2 to determine the local BHMF. All densities are computed 
with A(AfBH)=0.3. Uncertainties are the 16 and 84% percentiles of the Montecarlo realizations 
described in the text. These correspond to ilcr uncertainties. 



Adopted Luminosity A/bh^ Pbh (E+SO) pbh (All) 

or Velocity Function -host-gal. -prop. [ X 10^ Mq Mpc""^] 



Sheth et al. 2003 


Mbh-o-* 




5.0^1;! 


Bernardi et al. 2003 


A^BH--f'K,bul 


3-8t;:? 




Nakamura et al. 2003 


A/BH--^^K,bul 


2 5+^-2 
^•^-0.8 




Kochanek et al. 2001 


AfBH-^K.bul 


o n-l-1.0 
0.0_Q y 




Marzke et al. 1994 


A^BH-iK.bul 


3 4+2-7 


4.5+?:^ 



All L- or V- Func. ... ... 4.6 



tion o f the bulges of early types in the K band. ICole et all 
i200ll) estimate K — z* — 2.12 which, combined with 
the a verage value for the e arly type galaxies r* — z* = 
0.68 jB ernardi et al.l l2 nn3lJ). provides the required color. 
iFukugita. Shimasaku. fc Ichikawal (Il995l) estimate r* ~ z* = 
0.79 : 0.63 : 0.70 : 0.65 : 0.57 (E:SO:Sab:Sbc:Scd) and this 
justifies our conservative choice of the scatter which also 
allows us to apply the same color correction to all mor- 
phological types. For the A^BH-ibui relation we use the 
coefficients obtained by fitting the 'Group 1' galaxies of 
[Marconi & Hunt (2003) in the K band: 

logMBH = (8.21±0.07) + (1.13±0.12)(logLK,bui-10.9)(10) 

where LK,bui is in units of Lx©. The local BHMF derived 
with the MBH-ibui relation assuming intrinsic dispersions 
A(Mbh)=0.30 and A(Mbh)=0.5 are plotted in^. As be- 
fore, the 16 and 84% percentiles of 1000 Montecarlo realiza- 
tion of the local BHMF indicate ±la uncertainties (shaded 
area and errorbars). 

As clear from Fig. ^ the effect of including the in- 
trinsic dispersion in the AfBH-o"* and AfBH-^bui correlations 
is that of softening the high mass decrease of the BHMF, 
thus increasing the total density. But the most important 
result is that in order to provide the same BHMF, the re- 
lations Mbh-o"* and MeH-ibui mus t have the same intrin- 
sic dispersion to within 0.1 in log. IYu fc Tremaind ll2002l) 
found that pBH[AfeH — -^^bui] ~ 2pbh[Mbh — o"*] by using 
A(Mbh)=0 for Mbh-(t* and A(Mbh)=0.5 for MBH-Lbui- 
This discrepancy can be entirely ascribed to the efi'ect quan- 
tified by Eq. Indeed, the densities /9bh of the BHMF's 
plotted in ¥\g.^ are /9bh = 2.7l^;^ x lO'^ M© Mpc~^ [Mbh- 
ffi, wi th A(AfBH)=0, in agreement with IYu fc Tremaind 
(|2Qq3)] and pbh = 5.5li;3 x 10^ M© Mpc"^ [MBH-^bui 
with A(Mbh)=0.5]. Conversely, with A(Mbh)=0.30, one 
obtains iAt^l x 10^ M© Mpc"^ (A^bh-ct*) and 3.8li;g x 
l(f Mq Mpc"^ (A/bh-I/ bui) and these two values are in ex- 
cellent agreement. The densities in massive BH's were eval- 
uated in the log(AfBH/ M©) = 6—10 range and the same 
range will be considered throughout the rest of the paper. 

One advantage of using the A/bh-ct* relation is that 
the uncertainties on the derived BHMF are smaller than 
in the case of Mmi-Lhui- This is because with the AfBH-o-* 
relation one does not have to apply any correction for the 
bulge fraction. On the other hand, measuring stellar velocity 



dispersions is much more difficult than measuring galaxy 
luminosities; thus it is clear that the two relations A^bh-o"* 
and AfBH-i'bui should complement each other. 

In summary, the use of the same intrinsic dispersion 
for the A^BH-cr* and A/BH-ibui relations provides perfectly 
consistent BHMF's with the same mass densities pbh- This 
is a confirmation of the results bv iMarconi fc Hunt (200^ 
who showed that, when considering only secure BH measure- 
ments, the A^BH-cr* AfBH-ibui and A^BH-Afbui relations have 
similar intrinsic dispersions (~ 0.3 in log). When not taking 
into account the intrinsic dispersion of the Mbh-ct* rela- 
tion, the local BHMF is systematically underestimated at 
the high mass end, where the disagreement with the BHMF 
of AGN relics has been claimed. 

2.3 The Black Hole Mass Function for Early Type 
Galaxies 

Having established that both A/bh-o"* and AfBH-ibui rela- 
tions provide consistent BHMF's we can now evaluate the 
effects of using luminosity functions from different galaxy 
surveys and photometric bands in determining the BHMF 
in early type galaxies. 

In Fig. CId we compare the local BHMF's for early 
type galaxies obtained from different galaxy luminosity func- 
tions, in diffe rent photometric bands . We use the luminosit y 
functions bv iBernardi et al.l I 2003b^. iMarzke et al.1 |l994l . 
i Kochanek et alj J200ll) . and iNakamura et alJ ^|2003^ with 
details of the derivation specified in the following. 

• iBernardi et al.ll2003bl : this is the same BHMF plotted 
in p anel a derived using the A^BH-ibui relation. 

• iMarzke et al.lll994l : we use the luminosity functions per 
morphological type from the CfA survey. The luminosities 
are in Zwicky magnitudes, Mz, and we apply the color trans- 
formations directly measured by Eochanck ct al. ( 200j) 
computing Mz — K for all the objects used for the lumi- 
nosity functions {K is obtained from the 2MASS catalogue): 
Mz-K= (4.1±0.65;3.95±0.65;3.79±0.56;3.34±0.64) for 
(E;SO;Sa-Sb;Sc-Sd). The K magnitudes used are isophotal 
magnitudes 7^20 and the c orrection to total mag nitudes is 
Ktot = K20 - (0.2 ± 0.04) jKochanek et al.ll2001f) . Then we 
use the bulge-total correction and the AfBH-iK.bui relation 
fro m Tab. [J 

• iKochanek et alJl200ll : we use the luminosity function in 
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Figure 2. (a) Local BHMF's for all galaxy types derived using different galaxy surveys and shown with the same notation as in Fig.lTk. 
(b) Best estimate of the local BHMF obtained by combining all the Montecarlo realizations of the BHMF which were used for panel (a). 
The solid line represents the 50% percentile and the shaded area is delimited by the 16 and 84% percentile levels. 



the K band and we apply the correction from K20 to Ktot 
described previously, and we use the morphological fractions 
and Bulge- Total correction s from Tab.0 

• iNakamura et all 1200,4 we use the luminosity function 
in the SDSS r* band and we apply the color correction K — 
r* — (2.8 ± 0.2) described previously and the Bulge- Total 
corrections of Tab.Q 

As in the previous section, uncertainties are estimated with 
1000 Montecarlo realizations of the local BHMF where, in 
the case of the galaxy luminosity functions, we use only er- 
rors on $*, the galaxy number density. All the local BHMF's 
for early type galaxies are in remarkable agreement within 
the uncertainties. The discrepancy at the low mass end 
(Mbh < 10^ Mw) be tween the BHMF derived with the 
iBernardi et"all ^2003b^ luminosity functions and the others 
is not significant. It occurs in a region where the luminos- 
ity functions of early type galaxies are extrapolated and is 
due to the d ifferent functional form s adopted to fit the data 
(gaussian for lBernardi et al.l ll2003lJ) . Schechter functions for 
the others). 

2.4 The Black Hole Mass Function for All Galaxy 
Types 

Here we estimate the local BHMF by considering also late 
morphological types. We use all the luminosity and velocity 
functions describe d in the previous sectio n (with the excep- 
tion of the one bv IBernardi et a" ] J2003bll which is only for 
early type galaxies). For the velocity f unction of late type 
galaxies we take the estimate made bv ISheth et alj i2003fl 
and shown in their Fig. 6. When necessary, we apply the 
bulge-to-total correction with the numbers in Tab. Q and 
we then use either the MBH-^bui or the Mbh-o"* relation as 
described in Sec. 12.21 



Fig. shows the local BHMF for all galaxies while 
Tab. 121 reports the estimated BH mass densities pbh both 
for early type and all galaxies. These densities are com- 
puted for h = 0.7. All the BHMF's and /Obh's are in 
agreement within the errors. Finally, in Fig. |5|3 we present 
our best estimate of the local BHMF obtained by merg- 
ing all the random realizations of the BHMF's shown in 
Fig. and considering the 16, 50 and 84% percentile lev- 
els. Our best estimate of the local density in massive BH's 
is consequently pbh = 4.6lJ;4 ^ MqMpc"^. Roughly 
70% of the total BH density resides in early type galax- 
ies. Our e stimate of the density in loc al BH's is in agree- 
ment with iMerritt fc Ferraresd (1200181) and with lFerrares'3 
i2002h though, in the latter case, the shape of our BHMF is 
very different at the high mass end. O ur estimate is a fac- 
tor ~ 1.8 larger than th ose bv lYu fc^ rcmainc (2002) and 
lAller fc Richston3 (l2002ll and the reasons for this discrep- 
ancy are outlined in Sec.0] 



3 THE MASS FUNCTION OF AGN RELICS 

Once the mass function of local BH's has been determined, 
the subsequent step is the estimate of the BHMF of AGN 
relics, i.e. BH's in galactic nuclei which were grown exclu- 
sively during active phases from small (1-10^ Mq) seeds. We 
will first describe the continuity equation which will be used 
to relate the relic BHMF, N{M,t), to the AGN luminos- 
ity function (f){L,z), under the assumption that AGN are 
powered by mass accretion onto massive BH's. The physical 
quantity directly related to the mass accretion onto the BH 
is the total intrinsic AGN luminosity L. Since AGN luminos- 
ity functions are determined in limited energy bands, we will 
provide suitable bolometric corrections. We will then briefly 
describe the adopted AGN luminosity functions obtained 
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with our bolometric corrections. Finally we will present our 
estimates of the relic BHMF's derived from different AGN 
luminosity functions. 



3.1 The Continuity Equation 

The Black Hole Mass Function BHMF of AGN relics, 
N{M,t), can be estimated fr o m the continuity equation 
jCavaliere. Morrison, Woodl ISma,]l Bla!ndfo7dl 

Il992r with which, under simple assumptions, it is possible 
to relate N{M, i) to the AGN luminosity function. If 
N{M,t)dM is the comoving number density of BH's with 
mass in the range A'l and A'l + dM at cosmic time t, the 
continuity equation can be written as: 



dN{M, t) d 



dt 



+ g^[NiM,t){MiM,t))]=0 



(11) 



where {M{M,t)} is the "average" accretion rate on the 
BH of mass M. We adopt the working assumption that 
AGN's are powered by accretion onto BH's, and that the 
BH growth takes place during phases in which the AGN is 
accreting at a fraction A of the Eddington limit {L — ALecm) 
converting mass into energy with an efficiency e. We can 
thus simply relate the AGN Luminosity Function </>(!/, t) 
[0(1/, t)dlogL is the comoving number density of AGNs in 
the range log L, log L + d log L at cosmic time £] to the BH 
mass function; 



,^(L,t)dlogL = 5{M,t)N{M,t)dM 



(12) 



where S{M, t) is the fraction of BH's with mass M which are 
active at time t, i.e. the BH duty cycle. If a BH is accreting 
at a fraction A of the Eddington rate, its emitted luminosity 



(13) 



where Ie is the Eddington time, e is the accretion ef- 
ficiency and Mace is the matter falling onto the black 
hole. The growth rate of the BH, M, is thus given by 
M = {1 — e)Macc, since a fraction e of the accreted mass 
is converted into energy and thus escapes the BH. Since 
(M(M,t)) = 5{M,t)M{M,t), combining Eos. [HI and ITl we 
can write: 



N{M,t){M{M,t)} = 



(1 — e) r < , T M dL 

iq^[^(^'*)]-^^dM 

which can be placed in Eq. llll where the only unknown func- 
tion is N{AI,t). If e and A are constant we can then write 



dN{M,t) {l-e)\^c^ 



dt 



et% In 10 



d<t,iL,t) 



dL 



(15) 



which can be easily integrated given the AGN luminosity 
function and the initial conditions. 

For initial condition, we assume that, at the starting 
redshift Zs, S[M,t{zs)] = 1. This can be interpreted either 
by saying that at Zs all Black Holes are active or that we 
are following the evolution only of those BH's which were 
active at Zs- Thus 



MN{M,ts) = [0(L,t.)], 



(16) 



We will show that the final results are little sensitive to the 
choice of the initial conditions, provided that Zs > 3, since 
most of the BH growth takes place at lower redshifts. 

Eq. 1151 which represents the continuity equation in 
the case of constant e and A, can be trivially integrated 
on M to derive the relation used by various authors 
^Padflvanij Burg, fc Edelsonlll990l: IChokshi fc Tu7^ll992l: 
lYu fc Tremainell2002HFerraresjl2002li l 



PBH 



■Ut 



(17) 



note the factor (1 — e) which is needed to account for the part 
of the accreting matter which is radiated away during the 
accretion process. Ut is the total comoving energy density 
from AGN's (not to be confused with the total observed 
energy density) and is given by 



dz 



dt 
dz 



L(j)[L, z) dlog L 



(18) 



(j}{L, z) d log L is used if the AGN luminosity function is de- 
fined per logarithmic luminosity bin, otherwise it should be 
(j>{L, z) dL. 

The right-hand element of the continuity equation, 
which contains the source function, is null, meaning that 
we neglect any process which, at time t, might 'create' or 
'destroy' a BH with mass M. Indeed, in the merging pro- 
cess of two BH's, Ml -I- M2 ^ M12, BH's with Mi and M2 
are destroyed while a BH with M12 is created. We decided 
to neglect merging of BH's because, at present, the merg- 
ing rate, 'y{Mi,M2,z), is very uncertain and dependent on 
the assumptions of the model with which it is computed. 
Moreover, the aim of this paper is to assess if local BH's 
are indeed relics of AGN activity and by neglecting merging 
we can assess the importance o f mass accretion during AGN 
phases. IYu fc Tremaind (l2002l) have used an integral ver- 
sion of the continuity equation in order to take into account 
merging but their approach does not allow a direct compari- 
son of the local and relic BHMF's. Furthermore, the discrep- 
ancy they find between local BH's and AGN relics becomes 
worse if merging is important. iGranato et all i2004l) in their 
physical model of coevolution of BH and host galaxy find 
that the growt h of the BH i s due to mass accretion. Sim- 
ilarly, Haima n. Ciotti. fc Ostri kcr (2003i) find that merging 
of BH's is important only at high redshifts while, at lower z, 
BH growth is dominated by mass accretion. These reasons 
support our choice of neglecting the merging process. Indeed 
merging of BH's might well be very important in shaping the 
BHMF at high redshifts (2 > 3 — 5) but, as we will see, the 
relic BHMF at 2 = is very little dependent on its shape 
and normalization at 2 ~ 3. Moreover, if the relic BHMF at 
2 = is changed significantly by merging then its remark- 
able agreement with the local BHMF in the merging-free 
case would be a mysterious coincidence (see SecEJ. 



3.2 The Bolometric Corrections 

The BH mass accretion rate is given by the AGN luminosity 
L which, in turn, can be estimated from the luminosity in 
a given band b, Lb, by applying a suitable bolometric cor- 
rection fboi,b = L/Lb. At the beginning of Sec. |2lwe have 
referred to L as the total intrinsic AGN luminosity, and 
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Figure 3. (a) Template spectrum for a L = 10^^ Lq AGN obtained as described in the text (solid line). The dashed and dotted lines 
represent the radio quiet and radio loud templates by Elvis et al. 1994. (b) Bolometric corrections derived from the spectral templates 
built as described in the text. The error bars represent the 16 and 84% percentiles of the 1000 Montecarlo realizations described in the 
text. The hatched area represents the B band bolometric correction (with itlo" scatter) by Elvis et al. 1994. 



before presenting the adopted bolometric corrections it is 
necessary to define the meaning of intrinsic. 

The integral of the observed spectral energy distribution 
(SED) of an AGN provides the total observed luminosity, 
Lobs, thus the quasar bolometric corrections bv lElvis et alJ 
il994l) provide Lobs because they are based on the aver- 
age observed SED. However, Lobs does not give an accurate 
estimate of the BH mass accretion rate because it often in- 
cludes reprocessed radiation (i.e. radiation absorbed along 
other lines of sight and re-emitted isotropically) . The accre- 
tion rate is better related to the total luminosity directly 
produced by the accretion process, which we call the total 
intrinsic luminosity L. In AGN's, L is given by sum of the 
Optical-UV and X-ray luminosities radiated by the accre- 
tion disk and hot corona, respectively. Conversely, it is well 
kno wn that the IR radiation is reprocessed from the UV 
(e.g. lAntonuccilToQ^) . Thus, in order to estimate L one has 
to remove the IR bump from the observed SED's of unob- 
scured AGN's. In radio loud AGN's there might be an im- 
portant contribution from synchrotron radiation to the IR. 
However this will not affect our average bolometric correc- 
tions since the number of radio loud AGN's is small, ~ 10% 
of the whole AGN population. 

In previous works , several au t hors h ave used the bolo- 
metric corrections by lElvis et al.l (^O^) but, as explained 
above, these provide the observed and not the intrin- 
sic AGN lumin osity. Removing the IR bump from the 
lElvis et alJll994l spectral templates, the 1 /im-X-ray range 
includes ~ 2/3 of the total energy; thus one obtains, for 
instance, L/vbLvb = 7.9 ± 2.9 instead of the commonly 
used L/vbL^b = 11.8 ±4.3. However, the lElvis erallll994l 
quasar templates, obtained from an X-ray selected sam- 
ple of quasars, are X-ray "bright" (see the discussion in 



lElvis. Risaliti. fc Zamoranil |2002|) and thus underestimate 
the bolometric corrections in the X-ray energy bands. More- 
over the dependence on luminosity of the bolometric correc- 
tions cannot be obtained. For these reasons we will derive 
new estimates of the bolometric corrections as described in 
the following. 

First, we construct a template spectrum. In the optical- 
UV region our template consists of a broken power-law with 
Q = -0.44 [L^ (X u") in the range I < X < 1300 A 
and Q = -1.76 in the range 1200-500A (e.g. iTelfer et all 
120021: IVanden Berk et all 1200 j) . At A > 1 ^m our big blue 
bump is truncated assuming a = 2 as in the Rayleigh- 
Jeans tail of a blackbody. The X-ray spectrum beyond 1 
keV consists of a single power law plus a reflection compo- 
nent. Tlw_20werlawjias the typical photon index F = 1.9 
(e.g. iGeoree et aD Il998l: iPerola et all l2002l) and an expo - 
nential cutoff at Ec = 500 keV. FoUowing lUeda et al.H2003ll . 
the reflection component is ge nerated with the pexrav model 
iMagdziarz fc Zdziarskill995l) in the XSPEC package with 
a reflection solid angle of 27r, inclination angle of cos i = 0.5 
and solar abun dances. We then resca le the X-ray spectrum 
to a given aox iZamorani et alJll98J) which is defined as 

log[L42500A)/L42keV)] 

Ctox = 5^ (19) 

log[//(2500A)///(2keV)] 

and we finally connect with a simple powerlaw the point 
at 500A with that at 1 keV. To account for the de- 
pendence of Upx on lumino s ity, w e use the relation by 
IVignali. Brandt, fc Schneideil J2003D : 

Oox = -0.11 * log L^(2500A) + 1.85 (20) 

Finally, we assume that the template spectra, hence the 
bolometric corrections, are independent of redshift. 

The template spectrum for a. L — 10^^ L© AGN is plot- 
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ted in Fig.|^ and is compared with the rad i o quie t and radio 
loud quasar median spectra bv lElvis et al.l (Il994l) . Our tem- 
plate is in very good agreement in the optical-UV part but is 
obviously missing the IR bump, which is due to reprocessed 
UV radiation and thus not considered here. Our template 
is fainter in X-rays but this is expected si r ice. a s already 
discussed, the quasar sample bv lElvis et alJ il994l) is X-ray 
bright. 

In Fig.|Hj3 we plot the bolometric corrections L/vbLvb, 
L/L(0.5-2 keV) and L/L(2-10 keV) derived from the above 
templates. The errorbars lines represent the 16 and 84% 
percentiles from 1000 Montecarlo realization of the spectral 
templates where we have assumed the following la uncer- 
tai nties for the input parameters: ±0.1 for the s pectral slopes 
fa: lTelfer et al.ll2002l : IVanden Berk et al.ll200ll) in the 1 /^m- 
500A range and ±0.05 for the olox, constant at all lumi- 
nosities. The latt er is a conservative assumption (see, e.g., 
lYuan et al.lll99^ but is made to account for possible, but 
unaccounted for, systematic errors. As in the case of the lo- 
cal BHMF, the 16 and 84% percentiles correspond to ilcr 
uncertainties on the logarithm of the bolometric corrections. 
The 50% percentiles can be fit with a 3rd degree polynomial 
to obtain the following convenient relations: 

log(L/i(2-10 keV)) = 1.54 + 0.24/1 + 0.012/:^ - 0.0015/1^ 
log(L/L(0.5-2 keV)) = 1.65 + 0.22C + O.OUC^ - 0.0015-C^ 
\og{L/vBL^B) = 0.80 - 0.067-C + 0.017/1^ - 0.0023-C^ (21) 

where £ — (log L — 12) and L is the bolometric luminosity 
in units of L0. Hence, the bolometric corrections have log- 
normal distributions with average values given by the above 
equations and ±la" scatters (at fixed L) which can be derived 
from Fig.l^t- Scatters are given by ~ 0.05 for the B band and 
~ 0.1 for the X-rays, taken independent of L for simplicity. 
It is worth noting that th e B band bolomet ric correction is 
in agreement with that bv lElvis et alj (1199411 whose average 
value and scatter are shown by the h atched area in Fig.jSJj- 
The scatter of the lElvis et alJ 1^3) bolometric correction 
is the standard deviation of their quasar sample while ours 
are uncertainties on the average values. The average aox 
in the log(I// Lq) = 11.5 — 12.5 ran ge is 1.43, the same 
value estimated bv lElvis. Risaliti. fc Zamorani (2002) after 
correcting for the biases due to X-ray or optical selection of 
the parent quasar samples. 

3.3 The Luminosity Function of Active Galactic 
Nuclei 

In order to ensure a consistent treatment, the AGN lumi- 
nosity function 4>{L, t) used in the continuity equation must 
describe the evolution of the entire AGN population. The lu- 
minosity L is the total luminosity radiated from the accret- 
ing mass and is estimated with the bolometric corrections 
described in the previous section. 

The AGN luminosity functions available in the liter- 
ature do not describe the entire AGN population because 
they are the result of surveys performed in "narrow" spec- 
tral bands, with given flux limits an d selection criteria. For 
example, the luminosity function of iBovle et alJ ll2000h in- 
cludes only quasars selected for their blue color but it is 
now known there is a population of "red" quasars which 
is missed. The luminosity function of soft X-ray (0.5-2 



keV) selected AGNs of lMivaii. Hasinger. fc Schmidd (l2000h 
misses most of the sources with significant X-ray absorption 
{Nh > lO'^'^cm"^); it includes mostly broad lined AGNs 
(~ 80% of the total) but it is known that, at least lo- 
cally ^ there is a dominant p opulation of obscured AGN's 
(e.g. iMaiolino fc Riekelll995l) . Current hard X-rays surveys 
(2-10 keV) are less sensitive to obscuration than optical 
and soft X-ray ones , thus the very recent luminosity func- 
tion by lUeda et alj i2003h probes by far the largest frac- 
tion of the whole AGN population. However it is restricted 
only to Compton-thin AGNs {Nh < 10^* cm"^), while it is 
known that, at least locally, t here is a significant fraction 
of Co mpton-thick objects (e.g. iRisaliti. Maiolino. fc Salvatil 
^9^. In summary, when using the luminosity functions 
available in the literature, one should know that they de- 
scribe part of the AGN population, and thus account only 
for a fraction of the local BHMF. 

In this paper, w e consider the AGN luminosity functions 
bvl Bovle et all2000l fband -Bl.lMivaii. Hasinger. fc Schmidt 
2000^ (0.5-2 keV) and lUeda et al.ll2003l r2-10 keV). Thus, the 
AGN luminosity function which will be used in the continu- 
ity equation is given by 

(P{L) dL ^ (j>{Lx) dLx (22) 

where x is either the B, the 0.5-2 keV or the 2-10 keV band 
and L = fboi,x Lx- 

The AGN luminosity functions obtained with the 
bolometric corrections described in the previous sec- 
tion are compared in Fig. Hk, at select e d red shifts. 
The luminosity functions bv [ Bovle et al.l i2000f) and 
iMivaii. Hasinger. fc Schmidii (l200ol) are in rough agreement 
at the high L end, meaning that they may be sampling the 
same quasar population. Indeed most of the objects in the 
iMivai i. Hasinger. fc Schmidt (2000) sample are broad-lined 
AGN which , at h igh L, become the quasars observed by 
iBovle et al.l l|2000|). The disag reement at low luminosities is 
because the lBovle et al ] |2000l luminosity function is extrap- 
olated for A'Ib > —23 i .e. for L < 1 0^^ h p,. In contrast the 
luminosity function bv lUeda et alJ (|20o3) samples a larger 
fraction of the AGN population at all luminosities. 

In Fig. we compare the differential comoving en- 
ergy densities computed in the luminosity ranges 41.5 < 
log(Lx/erg§-i) < 48 (0.5 - 2keV and 2 - lO keV bands) 
and - 28 < Mb < -21 {B band). For the lUeda et al] 
i2003l) luminosity function we also plot the differential co- 
moving energy density for objects with L > lO^'^ Lq. 
High luminosity ob jects provide ^ 5 0% of the total en- 
ergy density in the lUeda et alJ (|2p03i) luminosity function 
and their emission has a redshift distribution similar to that 
of the AGN' s bvlMivaii. Hasinger. fc Schmidti ll2000tl and 
iBovle et al.l ll2000fl . Clearlv. lower luminositv objects con- 
tribute significantly at z < 1.5 and this is an imp ortant re- 
sult of the recent Ghandra and XMM sur veys fe.g. lHasingeJ 
bOOStlFiore et al.ll2003l : lUeda et alj|2003l) . The total comov- 
ing energy densities in the redshift range z = 0—3 (i.e. the in- 
tegrals of the quantities plotted in the fig ure) are 4 .2x 10~^^, 
6.4x10"^^ and l.SxlO"^ '^ ergcm~^ f or iBovle et al. 200'ol . 
iMivaii. Hasinger. fc Schm idt 20o3, a ndlUeda et al.l2003l! re- 
spectively. The high L objects in the llJed^^^dTl20o3) lu- 
minosity function provide 8xl0~^® ergcm"^, i.e. ~ 50% of 
the total energy density. 
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Figure 4. (a) AGN LF's at selected redshifts obtained with the bolometric corrections described in the text. Different symbols indicate 
different AGN luminosity functions, while different line styles distinguish between redshifts. (b) Differential comovi ng energy densities 
of AGN's from different surveys obtained with the bolometric corrections described in the text. In the case of the lUeda et al.l i2003h 
luminosity function we also plot the comoving energy density for high luminosity objects (L > 10^^ Lq). 'U', the total comoving energy 
density in the redshift range z = — 3, is in units of 10~^® ergcm"^ and '< z >' is the average redshift computed following Eg. 1291 




Figure 5. (a) Relic BHMF's at z = (solid lines) calculated assuming e = 0.1 and A = 1 compared with the relic BHM F's at Zs=Z 
(dotte d lines). Different symbols indicate the different AGN luminosity functions used, (b) Relic BHMF's obtained from the lUeda et alj 
i20n3r AGN luminosity functi on but with differe nt A, Zs and 5 with respect to panel (a). For comparison, the thick line with no symbols 
indicates the relic BHMF from lUeda et al .120031 plotted in (a). As in the previous panels, solid lines are the relic BHMF's at z = 0, while 
dotted lines are those at Zs- 



Local Black Holes, AGN Relics and the XRB 11 



By applying Eq. 1171 the mass densities in BH's can be 
written as 

PBH = p X 10^ Mo Mpc"^ (23) 

where p is 0.6 (from the iBovle et al.1 1200(1 AGN LF). 0.9 
jMivaii. Hasineer. fc Schmiddl200(i) ."2.2 JUeda et aLll2003ri 
and 1.2 (|Uedaet alJ 1 2003*. high L objects). The value 
derived from the lUed a ct al. (2003) luminosity function 
(2.2 X 10^ MoMpc"^) is already in marginal agreement 
with pBH = 4.5+i;4 X 10^ M© Mpc"^, i.e. the estimate of 
the local BH density given in Sec. 12.41 



3.4 The Relics of Active Galactic Nuclei 

Given the AGN luminosity func tions (|Boxle_et_alJ I2OO0I : 
iMivaii. Hasineer. fc Schmidtl2000HUeda et all2003ll we can 

integrate the continuity equation assuming that A = 1 
(AGN's emitting at the Eddington luminosity), e = 0.1 and 
2:5=3. Though we know that the above luminosity functions 
do not describe the whole AGN population, initially we do 
not apply any correction for the objects which have been 
missed. 

The relic BH mass functions are shown in F i g. 
As expected, the AGN's traced by the lUeda et alJ 11200311 
luminosity function leave more relics than t hose t raced 
by the luminosity functions of [ Bovle et alj ll2000h and 
iMivaii. Hasineer. fc Schmidt! ||200(J) . It is also notable that 
the relic BHMF's at z = (solid lines) are ~ 2 order of 
magnitude larger than those at = 3 (dotted lines) mean- 
ing that most of the BH growth took place for z < 2^ = 3. 
The dashed lines represent the relic BHMF obtained by con- 
sidering only AGN's with L > 10^^ Lq. [In practice this 
was obtained by multiplying the AGN luminosity functions 
by exp(-lOiVi)-] The comparison between the dashed and 
solid lines indicate that today high mass BH's (Mbh > 
10* M0) grew during quasar phases {L > 10^^ Lq). 

The BHMF's in Fig.|S^ were estimated assuming A = 1 
and e = 0.1. It can be seen from Eg. llSl that the efficiency e 
is a simple scaling factor and an increase in e will decrease 
the level of the relic BHMF. However a variation of the Ed- 
dington fraction A will not have a simple scaling effect but it 
will produce a combination of scaling and translation along 
the M axis since A also enters the L — XM c? /Ie relation. In 
Fig.Et" the same relic BHMF derived using fUeda et al]l2003l 
(thick line with no symbols) is compared with the BHMF 
obtained assuming A = 0.1 (line with empty squares). As 
previously, solid and dotted lines indicate the BHMF at 
2 = and at 2 = 2s respectively. Decreasing A has the 
net effect of increasing the number of BH's at high masses 
(M > 10®'^ M0) while decreasing that at lower masses. 
In the same figure, we also show the BHMF's obtained by 
assuming that the starting redshift of integration is 2s = 5 
(line with empty triangles) and 2s = 2.5 (line with empty di- 
amonds) and that the fraction of active BH's at 20 is 5 = 0.1 
instead of 1 (line with stars). The relic BHMF at 2 = is 
almost independent of the initial conditions. Changing the 
starting assumption (zs, or 5[zs\) does not have any appre- 
ciable effect provided that 2s > 2.5. This is a consequence 
of the fact that the main BH growth takes place at 2 < 3 
when there is more time available (~ 85% of the age of the 



universe). At larger redshifts there is too little time for the 
BH growth. 

In summary, relic BH's grew mainly at low redshifts 
(2 < 3) and high mass BH's were produced during quasar 
activity. The initial conditions for the continuity equation 
influence little the relic BHMF at 2 = 0. 



4 LOCAL BH'S AND AGN RELICS 

In this section we compare the local and relic BHMF's. At 
first we will consider the relic BHMF's derived in the pre- 
vious section, without trying to account for the AGN pop- 
ulation missing in the adopte d luminosity func tions. After 
showing that, as expected, the lUeda et al.l l|2003h luminosity 
function encompasses the largest fraction of the AGN pop- 
ulation, thus providing a better match to the local BHMF, 
we will focus on it and try to account for the missing AGN's 
in a way which also satisfies the constraints imposed by the 
X-ray background. 

In Fig. Hi we compare the local BHMF derived in 
Sec IT^ (Fi g. E h) with the relic BHMF's obtained adoptiii g 
the lBwle etIS] j2000MMivaii. Hasineer. fc Schmidl J200CI) , 
and lUeda et alJ J2003D AGN luminosity functions (e = 0.1, 
A = 1) without any correction for the missing AGN popula- 
tion. As shown in the previous section these relic BHMF's 
are insensitive to the adopted initial conditions thus the pa- 
rameters they depend on are only e and A. 

The first result which can be evinced from the figure 
is that there is no discrepancy at high masses. All the relic 
BHMF's are slightly smaller than the local BHMF. As ex- 
plained previously, this is not significant because all the 
adopted luminosity fu nctions are missing pa rt of the AGN 
population. Conve rselvlYu fc Tremaind ll2002l) and lFerraresd 
(2002) using the lBov'l^^an ~ il200ol) luminositv function 
found a significant discrepancy because their relic BHMF 
at high masses is larger than the local one. The reasons for 
the solution of this discrepancy are the following. 

• We have taken into account the intrinsic dispersion of 
the AfeH-o"* and MBH-^bui correlations thus softening the 
high mass decrease of the local BHMF. 

• We have adopted zero points in the Mbh-o"* and A/bh- 
Lbui correlations which are a factor ~ 1.6 ( ~ 0.2 in log) 
larger than those used by lYu fc Tremaind ll2002ll . These 
zero points (and slopes) were derived from the analysis of 
^arconi & Hunt (2003) and their larger value is a con- 
sequence of the rejection of galaxies with unreliable BH 
masses. 

• We have adopted bolometric corrections which, on av- 
erage, are ~ 2/3 of those adopted by previous authors. This 
is because we have not taken into account the IR emission 
in the estimate of the bolometric luminosity. The IR bump 
is produced by reprocessed UV radiation and thus its use in 
the determination of L would result in overestimated accre- 
tion rates. 

We now focus on the lUeda et alJ (1200311 AGN luminosity 
function since we have established that, as expected, it is 
the one which provides the most complete description of 
the AGN population. The aim is to assess whether the relic 
BHMF can account for the whole local BHMF, i.e. if local 
BH's can be entirely explained as being AGN relics. 
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Figure 6. (a) Local and relic BHM F's. The shaded ar ea represents the itlcr uncertainty on the logarithm of local BHMF (see Fig. 1^). 
(b) Relic BHMF obtained using the lUeda et al.l l2003^ hard X-ray luminosity function and accounting for the missing Compton Thick 
AGN's (see text). The errorbars are the 16 and 84% percentile levels of 1000 realizations of the relic BHMF (see text). The local BHMF 
is plotted with the same notation as in the panel (a). 



The lUeda et alj ll2003h luminosity function only de- 
scribes the population of Compton-thin AGN's, i.e. ob- 
jects with absorbing column densities logA''£r < 24 [cm~^]. 
To estimate the correction for the missing Compton-thick 
sources, we follow lUeda et al.1 (|20o3) and assume that there 
are as many AGN's in the 24-25 bin of log Nh as in 
the 23-24 one, as indicated by the Nh distribution of 
iRisaliti. Maiolino. fc Salvatilll999l Thus, w ith the Nh dis- 
tribution estimated bv luede^t al.l ll2003l) . the luminosity 
function has to be multiplied by 1.3 in order to account 
for the Compton-thick AGN's. In reality, the correcting 
factor is luminosity dependent but, since it is limited to 
the 1.2 — 1.4 range, we chose an average value for sim- 
plicity. W ith this contribution from Compton-thick sources, 
lUeda et' al. (2003) are able to provide a reasonable fit of the 
X-ray background spectrum. (We will return to this issue 
in the following section.) While AGN's at log Nh > 25 are 
not important contributors of the X-ray background, they 
do contribute to the relic BHMF. Thus, we assume that 
for log Nh > 25 there are as many AGN's as in the 23- 
24 or 24-25 log Nh bins. This assumptio n is also justified 
by the lRisaliti. Maiolino. fc Salvatil il999l) Nh distribution. 
The correcting factor becomes then 1.6. Therefore, to com- 
pensate for the missing obscured AGN's we apply a correct- 
ing factor of 1.6 independen tly of luminosity. In Fig. |Sj3 the 
corrected lUeda et alJ i2003h luminosity function is used to 
determine the relic BHMF. We also show the uncertainties 
(16 and 84% percentiles, i.e. ±1(t errors on the log of the relic 
BHMF as previously) estimated with the usual 1000 Mon- 
tecarlo realizations of the relic BHMF. These were obtained 
by varying the number density of the luminosity function 
(±30% la error on the number density to avoid correlated 
errors on the other parameters), the hard X-ray bolometric 



correction (±0.1 la error on log fboi,x), and the factor used 
to correct for the missing Compton-thick AGN's (±0.3 la 
error). 

The relic and local BHMF agree well within the uncer- 
tainties. Thus, adopting the best possible description of the 
whole AGN population, the mass function of relic BH's is 
in excellent agreement with the mass function of local BH's. 
Local BH's are thus AGN relics and were mainly grown dur- 
ing active phases in the life of the host galaxy. This agree- 
ment has been obtained with e = 0.1 and A = 1 indicat- 
ing that (i) efficiencies higher than commonly adopted for 
AGN's are not required, and (ii) the main growth of BH's 
occurs in phases during which the AGN is emitting close to 
the Eddington limit. In Section|S]we will explore the locus in 
the A — e plane which is permitted by the comparison from 
the local and relic BHMF's. As discussed in Section f3. II we 
have neglected merging in our estimate of the relic BHMF. 
However, the good agreement between the local and relic 
BHMF's suggests that the merging process does not signifi- 
cantly affect the build-up of the BHMF, at least in the z < 3 
redshift range. 



5 CONSTRAINTS FROM THE X-RAY 
BACKGROUND 

From the X-ray background (XRB) fight it is pos- 
sible to estimate th e expe ct ed mass density of relic 
BH's ("Salucci et al.' '1999'; ' Fabian fc Iwasawal 119991 : 
^Elvis. Risaliti^ fc Za m orani 2002^) which can then be 
compared wit h the mass density of loca l BH's. From this 
comparison, lElvis. Risaliti. fc Zamoranil ll2002t) inferred 
that massive BH's must be rapidly rotating for the high 
efficiency needed (e > 0.15) to match the XRB and local 
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BH mass densities. We will show that, due to the redshift 
distribution inferred from the new X-ray surveys, the match 
between the XRB and local BH mass densities can be 
obt ained without requiring large efficiencies (e > 0.1; see 
also iFabianlEool IComastril 12003) . 

The XRB provides anot her type of con- 
straint. It has been shown (|Sefrti^Jg_Woltiej ^83; 
Madau. Ghisollini. & Fabian '19941: IComastri et alJ 1199,4 
Gilli. Risaliti. fc Salvati .1999.1 that the XRB spectrum can 



be reproduced by summing the spectra of the whole AGN 
population after suitable corrections to take into account 
the absorption along the line of sight. The AGN population 
and its redshift distribution is derived from AGN luminosity 
functions but, in order to fit the XRB spectrum, one has to 
include a correction for the missing (obscured) AGN's. This 
correction is the same one which must be adopted here to 
determine the relic BHMF from the whole AGN population. 
Thus, the ratio R\2 between AGN's included and missed 
in the adopted luminosity functions is a parameter which 
determines both the relic BIfMF and the XRB spectrum 
(and the X-ray source counts). When rescaling the AGN 
luminosity function with the (1 -f Ru) factor to match the 
local BHMF, one should also verify that it is also possible 
to match the XRB spectrum and source counts at the same 
time. In practice, one can use the R\2 value required by the 
XRB synthesis models in order to fit the XRB spectrum 
and source counts. It is beyond the scope of this paper to 
produce a model synthesis of the XRB but we will show 
that with 7?i2 values found in the literature we can match 
the local and relic BHMF, and also satisfy the constraints 
imposed by the X-ray background. 

The local density in massive BH's expected from the 
observed X-ray background (XRB) light can be estimated 
with the relation 



1 



PBH 



(1 + {z)) m 



(24) 



where is the total observed (as opposed to comoving) 
AGN energy density and {z) is the average source redshift. 
The factor (1 — e) is needed to take into account the fact 
that not all the accreting mass falls into the BH. C/y can be 
estimated from the observed X-ray background light as 



* 47r 

Ut ~ fbol.X X fobsc X Ix 

C 



(25) 



where Ix is the total observed surface brightness of 
the X-ray background (i.e. the integral of the XRB 
spectrum), fotsc is the correction to take into ac- 
count source obscuration in the X-rays (i.e. fobsc x Ix 
would be the total XRB surface brightness if AGN's 
were not obscured) and fboi,x is th e X-ray bolomet- 
ric correction (see for more details ^^im^^^wasawal 
I1999I; 'Salucci et al."l999'; 'Elvis. Risal iti. fc Zamoranill2002l) . 
|Elvis. Risaliti. & Zamorani (2002) estimate ~ (1.5 — 
3.4) X lO^^^ergcm"^. 

We first verify that the above formula is consistent with 
the scheme followed in this paper, and then establish how 
(z) must be computed. The observed background surface 
brightness at energy E is 



4nDl 



dz 



f[E{l + z)]L(l>{L,z)dL (26) 



where Dl is the luminosity distance, V is the comoving vol- 
ume, L is the source luminosity in the energy band x and 
(j) is the luminosity function in the same band. f[E] is the 
source spectrum at energy E normalized to have unit lu- 
minosity in the band x. Integrating on E to find the total 
surface brightness in band x one finds 
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(27) 



Applying the bolometric correction [L — fboi.xL with 
(j}{L,z)dL — <j){L,z)dL), the obscuration correction (fobsc) 
and comparing with Eq. 1181 one finds that 
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(1 + {z))U} = Ut 



(28) 



where It = fboi,xfobscIx and Ut is the total comoving (as 
opposed to observed) energy density. The average redshift 
(z) is then 
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' (l + z) dz 

Using the luminosity function by lUeda et alJ |2(i 
one finds (z) — 1.1, lower than the value (z) ~ 2 
assumed by previous authors. For compariso n, the high 
lumin osity objects (L > 10^^ L0) in the lUeda et alJ 
1200^ luminosity function have (z) = 1.3 while the 
iMivaii. Hasinger. fc Schmidt (20oJ) and lBovle et al. (200^ 
have (z) = 1.4 and (z) = 1.5, respect i vely. With the Ut 
estimate bv lElvis. Risaliti. fc Zamoranil ll2002t) we get 



PBH = (4.7 - 10.6) 



(1- 



X 10^ M0 Mpc" 



(30) 



which is perfectly consistent with the estimate from the lo- 
cal BHMF pBH = (3.2 - 6.5) x lO'^ Mq Mpc"^ without re- 
quiring efficiencies larger than the 'canon ical' value e= 0.1. 
This agreement has also been remarked bv lFabi;inll2003l and 
lGomastrilf2003l . The critical point is clearly the value of (z), 
the average redshift of X-ray sources emitting the XRB, 
which has be en sign ificantly r educed by Chandr a and XMM 
surveys Ce.g. iHasiuf^ er 2003: .Fiore et al. 2003,: Steffen et alJ 
l2003l:IUeda et aU2003D . The minimum value of e, allowed for 
a consistency between the two estimates of Pbh, is e = 0.07 
just slightly larger than the non-rotating BH case. Con- 
versely, the maximum allowed efficiency is e = 0.27, sub- 
stantially below the maximally rotating Kerr BH case (see 
Sec. |SJ . It is intriguing to find that efficiencies smaller than 
those expected from non-rotating BH's or larger that those 
expected from maximally rotating Kerr BH's are excluded. 

We have thus verified that the expected density of BH 
remnants inferred from the XRB is consistent with the lo- 
cal one without requiring efficiencies larger than the canon- 
ically adopted value e = 0.1. We now verify if, with the 
obscured/unobscured ratios R12 adopted in XRB synthesis 
models, it is possible to reproduce also the local BHMF. 
For the scope of this paper it suffices to notice that the 
lUeda et al.l l|2003r) luminosity function, with the correc- 
tion for the missing Compton-thick AGN's that we also 
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adopt, is used by the same authors to successfuUy repro- 
duce the XRB spectrum and source counts. Thus the agree- 
ment of the local and the relic BHMF in Fig. is ob- 
tained by also meeting the constraints from the XRB which, 
in practice, provide an estimate of the number of AGN's 
missed by the luminosity fun ction. The same comparison 
could also be done using the iMivaii. Hasinger. fc Schmidti 
\20Q ^ lumin osity function combined with the background 
model of Gill i. Salvati. fc Ha singer ( 20Q1.) . However, though 
that model is successful in reproducing the XRB spectrum 
and source counts, recent Chandra and XMM surveys have 
sho wn that the mo del redshift distribution is not correct 
fe.g. lHasingeill2003l) . 

In summary, our analysis shows that it is possible to 
meet the XRB constraints both in terms of pbh and of the 
local BHMF. 



ACCRETION EFFICIENCY AND 
EDDINGTON RATIO 



The rehc BHMF derived from the lUeda et alj (l2003l) 

AGN luminosity function, corrected for the Compton-thick 
AGN's, provides a good match to the local BHMF and also 
satisfies the XRB constraints. The match is obtained for 
e = 0.1 and A = 1. Here we investigate the locus in the 
A — e plane where an acceptable match of the local and relic 
BHMF's can be found, i.e. we determine the acceptable A 
and e values. 

To quantify the comparison between the local and relic 
BHMF's we recall that all the reahzations of the BHMF's 
have a log-normal distribution at given M and we consider 
the following expression: 



fc2= / LlogiVLW-logiVfl(M)]- ^^^^^^^ 



dlogM(31) 



where Nl{M) and Nr{M) are the local and relic BHMF's 
with their Icr uncertainties, ul^M) and aFi{M). The inte- 
gration is performed in the log(A4"/ Mq) =6—10 range, 
fc^ is the average square deviation between the logarithms 
of the two BHMF's measured in units of the total stan- 
dard deviation. If, for instance, logA'^L(Af) = log Nii{M) + 
ny''aL{M)'^ + aii{M)^, then — ri^, i.e. the two functions 
differ, on average, by n times the total standard deviation. 
The e and A values that corresponds to the minimum 
(^min — 0.4^) are marked by the filled square and are 
e = 0.08 and A = 0.5. To have an acceptable match be- 
tween the local and relic BHMF we require that ^ 1 
and this constraint identifies the region limited by the solid 
line in Fig. |7| For comparison, the dashed line limits the 
region where k^ < 0.7^, which corresponds to the aver- 
age square deviation in the "canonical" case, e — 0.1 and 
A = 1, which is marked by the cross. The allowed efficien- 
cies include the non-rotating Schwarzschild BH (e = 0.054, 
ISchwarzschildllQldlShapiro fc Teukolskvll983D and are well 
below the maximally rotating K err BH (e = 0.42, iKerJ 
Il963l: IShapiro fc Teukolskvlll983l) . The best agreement be- 
tween the local and relic BHMF's is obtained for efficiencies 
larger than that of the non-rotating B H case suggesting that, 
on av erage, BH's should be rotating. iHughes fc BlandfordI 
i2003h found that BH's are typically spun down by merg- 
ers and this would limit the importance of mergers in the 
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Figure 7. Locus where accretion efficiency e and Eddington ratio 
A provide the best match between the local and relic BHMF's. 
The solid and dashed lines limit the regions where the logarithms 
of the local and relic BHMF difler, on average, by less than 1 
and 0.7a (see text). The filled diamond marks the minimum 
and corresponds to e = 0.08 and A = 0.5. The cross indicates 
e = 0.1 and A = 1. The dashed lines are theoretical efficiencies for 
non-rotating (e = 0.054) and maximally rotating BH's (e = 0.42). 



growth of BH's, in agreement with our assumption in the 
continuity equation. 

The allowed Eddington ratios, A = L/L-Edd, are in the 
range 0.1 < A < 1.7 indicating that BH growth takes place 
durin g luminous accretion phas es close to the Eddington 
limit. iMcLure fc Dunlotil ll2003l) . using a large sample of 
SDSS quasars, have recently estimated that the average A 
varies from 0.1 at z ~ 0.2 to 0.4 at z ~ 2. Accounting for 
the different bolometric corrections (they used fboi.B ~ 9.8, 
~ 1.5 times larger than the value adopted by us), A varies 
from 0.15 at z ~ 0.2 to 0.6 at 2; ~ 2. These values are in ex- 
cellent agreement with the constraints posed on A in Fig. |7| 

The results in Fig. |7|do not imply that accreting BH's 
cannot have e and A values outside the region of the best 
match between the local and relic BHMF's. Indeed, those 
limits are only for average efficiencies and Eddington ratios 
during phases in which Black Holes are significantly grown. 

We conclude by noting that e represents only the ra- 
diative efficiency. If there is significant release of mechan- 
ical energy, the true efficiency might be higher. For in- 
stance, in M87 the kinetic energy c arried away by the jet is 
much larger than the radiated one iOwen. Eilek. fc KassimI 
I2OOOI) and, in general, the jets of radio loud AGN's can 
carry away up to half of the total power in kinetic energy 
iRawlings fc Saunders*199l':! Celotti. Padovani. fc Ghisellinil 
1I997,; Tavecchio et al. _200Q(). The release of mechanical 
and radiative energy is important for the feedback on the 
galaxy which is thought to be one of the causes behind the 
Mbh-ctj, and MBH-J^bui correlations (e.g. |Silk fc Reejll998l : 
lBlandfordl[l999l : lBegelmanll2003l:lGranato et alJl20ol) . Tak- 
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Figure 8. (a) Growth history of total BH mass density (with e = 0.1 and A = 1). Different symbols indicate the use of different AGN 
luminosity function s. The dashed and dot-dashed lines indicate the growth history of pbh in low < 10* Mq) and high mass BH's, 
computed using the lUeda et al.l 1200,^) luminosity function. The densities are normalized to the final value at 2: = which is reported in 
the upper part of th e plot. The units f or pbh 10^ Mq Mpc~^. (b) Average growth history of BH's with given starting mass at Zg=3 
computed using the lUeda et al.l 1200.^) luminosity function and e = 0.1, A = 1. The symbols, filled squares, empty squares , and stars, 
indicate when the BH reaches 90%, 50%, and 5% of its final mass, respectively. 



ing into account mechanical energy, the expression — e) 
should be transformed to — en — em), where en and 

Em are the radiative and mechanical efflciencies. Then, sim- 
ilarly, to what has been done in section, one could place 
constraints on both tR and em but this is beyond the scope 
of this paper. 



7 GROWTH AND ACCRETION HISTORY OF 
MASSIVE BLACK HOLES 

Having established th at the obscuration corrected AGN lu- 
minosity function by lUeda et al.l ^200.^1 with e = 0.1 and 
A = 1 provides a relic BHMF which is fully consistent with 
the local BHMF and the X-ray background, we can analyze 
the growth history of massive BH's and, in the next section, 
the average lifetime of their active phases. 

The redshift dependence of the total density in massive 
BH's is given by CEas. ITTland lll^ 



PBH (2') 



dz 



dt 
dz 



L(f)(L, z)d\og L 



(32) 



and is plotted in Fig. The rel ics of the AGN's t raced 
with the luminosity functions bv iBovle et al.l ll2000ll and 
iMivaii. Hasinger. fc Schmidt! (|200nl ) reach 50% of the 2 = 
mas s density around 2 ~ 1.7, while the AGN's traced by 
the lUeda et al.l ll2003l) luminosity function do so at 2 ~ 
1.4. This is a consequence of the larger number of low 
luminosity AGN's which are present at 2 ~ 1 (see also 
Fig. Bt*)- Indeed, when separating the contributions from 
low (A/bh < 10* M0) and high mass BH's which contribute 
50% each of the 2 = BH density, it is clear that low 



mass BH's grow later than high mass BH's. T hese, like in 
the case of the AGN's traced bv iBovle et alj (igOOQ) and 
iMivaii. Hasinger. fc Schmidd J2000l) . reach 50% of their fi- 
nal mass at 2 ~ 1.7. 

This issue can be further investigated by computing the 
average growth history of a BH with given starting, or final, 
mass. The growth history of a BH with given starting mass 
Mbho at Zs can be estimated as follows. The average accre- 
tion rate at 2 (or t) is given from Eq. 15 as: 



{M(M,t)) = 



(l-e)A 



ffi InlO eN[M,t) 



(33) 



thus one can solve the following differential equation to ob- 
tain the 'average' growth history of BH's: 



dM = {M(M, t)) dz 
dz 



(34) 



In Fig. IHJj we plot the average growth history of BH's with 
different sup ermassiv e BH like that 

of M87 or Cygnus A (Mbh ~ 3 x 10^ Mq; iMarconi et"al] 
I1997I: iTadhunter et aljT2003D was already quite massive (~ 
10^ M(7)) at 2s wh ile a smaller BH like that of Centaurus 
A iMarconi et akl lioOl ) was less massive, around A/bh ~ 
10^ Mq. a supermassive BH with Mbh ~ 10® Mq at 2^ = 3 
should now be over 10^° Mq. Indeed, the existence of very 
massive BH's at high redshifts is suggested by the detection 
of very luminous quasars and, in particular, those detected 
at 2 ~ 6 by the SDSS survey (e.g. Fan 2003). For instance, 
A^BH for the farther quasar known (2 ~ 6) is estimat ed as 
Mbh = 3 X lO'^ Mq dWillott. McLure. fc JarvT^bOOSH . and 
one would expect its local counterpart to be more massive 
than 10^^° Mq. However, these quasars have not been de- 
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Figure 9. Cosmic BH accre tion history compared with the esti- 
mate by|Barger_et_alJ ^OOll) and with the cosmic star formation 
rate bv ICharv fc EIbad^2001^ . 



tected yet and it is not clear if this is because these hyper- 
massive BH's are very rare or they simply do not exist. In the 
latter case there should be a physical reason which prevents 
a BH from growing beyond 10^" M©, possibly the feedback 
on t he host galax y mentioned in the previous sections (see 
also lNetzeil l2003f). 

From Fig. EJa we can also infer a confirmation of what 
already found in Fig. namely that more massive BIf's 
grow earlier. The symbols in the figure, filled squares, empty 
squares, and stars, mark the points when a BH reaches 90%, 
50%, and 5% of its z = mass, respectively. It is clear that 
for z < Zs = Z a\\ BH's gain more than 95% of their final 
mass but BH's which are now more massive than 10* M© 
had already gained 50% of their 2 = mass at z ~ 2. Con- 
versely, BH's which have now masses around A^bh ~ 10^ M© 
grew very recently, at z < 1. Again, this i s a co nsequence 
of the luminosity function by iUeda et al.1 J2003ft in which 
the distribution of lower luminosi ty AGN's p e aks at 2; ~ 1. 
Thus, the luminosity function of lUeda et al.l (|200^ points 
toward an anti-hierarchical growth of BH's in the sense that 
the largest BH's were formed earlier. 

If the correlations between Mbh and host galaxy prop- 
erties were v alid at higher reds hifts (as is suggested by 
the results of IShields et alJ |2003^ this would immediately 
imply that also the most massive galaxies should form 
earlier, in contrast with the predictions of current semi- 
an alytic models of galax y formation (see the introduction 
of iGranato et alj i2004l) for more details and references) . 
The detection of high mass galaxies in sub-mm surveys is 
indeed more consistent with the 'monolithic' scenario in 
whic h massive elliptica l s form at relatively high redshifts 
(e.g. iGenzel et all 120031: iGranato et aPbOO^I and references 
therein) . 

The BH accretion history (i.e. the total accretion rate 



at given z per unit comoving volume) can be estimated using 
Eg. in as 

/•M2 

Magn{z)= N[M,t{z)]{M[M,t{z)])dM (35) 

J Ml 

In Fig. |5| we plot the cosmic accretion history of BH's 
(computed using the lUeda et al.l (12003) luminosity function, 
e = 0.1 and A = 1) and we compare it with the estimate by 
|Barge^^alJjj20^)) and with the cosmic star formation rate 
bv lCharv fc ElbazI (20011. 

Our estimate of th e cosmic accretion rate onto BH's 
agrees well with that bv lBarger et alJ (120011) . likely because 
most of the hi j^h redshift AGN's u sed bv lUeda etall (l2003l) 
come from the'Barger et alJ (1200 sample. Differently from 
us,|BarEcr et al. (2001) estimate the bolometric AGN lumi- 
nosities by integrating the observed spectral energy distri- 
butions. Thus the agreement with our analysis should be 
viewed as a consistency check on the bolometric corrections 
and on the corrections for obscured AGN's that we applied. 

The cosmic accretion history has a similar redshift de- 
pendence as the cosmic star f ormation rate, which w e report 
here in the form estimated bv lCharv fc ElbazI j200lD . To aid 
the eye, the dashed line represents the estimate of the cosmic 
accretion history rescaled by 4000. The comparison suggests 
that indeed the two rates have a similar redshift dependence 
and justifies the assumption that the accretion onto a BH 
is proportional to the star formation rate, at least at a cos- 
mic level. This fundamental assumption is made in several 
models of coeval evolution of BH and galaxy and, together, 
with the feedback from the AGN expl ains the observed cor- 
relations Mbh-ct-i, and MBH--t/bui (e.g. IGranato et aljl20o3 : 
iHaiman. Ciotti. fc Ostrikerll2003fl . 



8 THE LIFETIME OF ACTIVE BH'S 

We now estimate the average lifetime of active BH's with 
the formalism used in this paper and t he corrected AGN 
luminosity function bv lUeda et alJ l|200^. 

In Fig. we plot the BH duty cycle 5{M, t) at se- 
lected redshifts de rived from Ea. ll2l and computed using the 
lUeda et al .* f2003') luminosity function, e = 0.1 and A = 1. 
From Sec. 13. II the duty cycle 5{M, t) is the fraction of BH's 
with mass M active at time t or redshift z. According to 
the definitions used in this paper, we consider a BH active 
if it is emitting at the adopted fraction A of the Eddington 
luminosity. Hence, objects which are usually classified as 'ac- 
tive' but which are emitting well below their Eddington limit 
(e.g. M87 or Centaurus A), should not be counted among 
the active BH's whose fraction is given by the duty cycle. 
BH's more massive than ~ 10^ Mq are very rarely active 
in the local universe (only 1 out of 10000) while they be- 
come more numerous at higher redshifts (by a factor ~ 100 
at 2; = 2). Conversely, lower mass BH's are usually a factor 
10 more numerous. A unit duty cycle at z = 3 is the initial 
condition assumed for the solution of the continuity equa- 
tion (see Sees. and The values of the duty cycle 
we obtain at z ~ — 1 are cons istent with the average val- 
ues of 3 — 6 X 10~^ estimated bv lHaiman. Ciotti. fc Ostrikerl 
J2003h . 

The average duration of the accretion process, i.e. the 
mean lifetime of an AGN which has left a relic of mass M, 
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is estimated by first solving Eg. 1341 to obtain M{z, Mo), the 
growth history of a BH with mass Mq at Zs- Then the 'active' 
time is simply given by 

, dt 



t{M) 



S[M{z,Mo) 



dz 



-dz 



(36) 



In Fig. If Ob we plot the average mean lifetime of AGN's (solid 
line) as a function of the relic BH mass &t z = computed 
using the lUeda et alJ |2003i ) AGN luminosity function and 
e = 0.1, A = 1. The average lifetime of AGN's which leave a 
relic BH mass of > 10^ Mq is of the order of 1.5xl0*yr, 
while for smaller relic masses (< 10* Mq) longer active 
phases are needed {tbh ~ 4.5 x 10* yr). Considering two lim- 
iting cases from Fig. □ (e, A) = (0.04,0.1), and (0.15,1.6), 
the average lifetimes can increase up to lO" yr. We remark 
that these numbers are the average lifetimes for z < Zs with 
Z3=3. 

The derived lifetimes are consistent with the main re- 
sult from this paper that most of the BH masses are assem- 
bled via mass accretion. Indeed, if a BH grows by accreting 
matter with efflciency e and emitting at a fraction A of its 
Eddington lu minosity, its e-f olding time is given by ts, the 
Salpeter time JSalpeteilll963) . 



(l-e)A 



= 4.2 X 10^ 



9e 



A" 



(37) 



Thus, the initial BH mass has been e-folded ~ 3 times from 
Zs = ?< for A/bh > 10^ Mq and more than 10 times for 
Mbh < 10* Mq. The apparently long lifetimes are thus the 
natural consequence of the fact that to grow a BH from small 
(~ 1-10 Mq) or iriterme diate mass seeds (100 - 1000 Mq, 
e.g. ISchneider et al.ll2002l) several e-folding times must pass 
easily implying tbh > 10* yr. For instance, to grow a BH 
from 10^ Mq to 10^ Mq with A = 1 and e = 0.1, one would 
need ~ lits i.e. ~ 6 x 10 * vr (see also|H a man fc Loebl200il) . 
Indeed, in the model bv lGranato et alf i 2004^ of coeval evo- 
lution of BH and galaxy, the time needed to grow most of 
the BH mass is ~ 3 x 10* yr (note however that they use 
e = 0.15, A = 3 with which the Salpeter time is roughly 
halved with respect to our paper). 

The reason why the growth of smaller BH's (< 10* Mq) 
requires longer time is not obvious but can be understood 
from Fig. Consider, for instance, BH's with masses ~ 
10^ Mq and ~ 10^° Mq at z = 0. Their average growth 
history is traced by the lower and upper tracks in the figure. 
The larger BH has Mbh ~ 10^ Mq at z = 3 and has to 
e-fold its mass by 3 times. Conversely, the smaller BH has 
Mbh ~ 10^ Mq at z — 3 and has to e-fold its mass by at least 
ten times. The masses of large BH's (> 10* Mq) increase 
by smaller factors from their z = 3 values, compared to the 
masses of smaller BH's. Thus they require shorter active 
phases. 

In general, literature estimates of AGN lifetimes range 
from 10 ^ to 10* yr and are still much uncertain (see the re- 
view bv lMartinill2003r) . Models where BH's grow by a com- 
bination of gas accretion traced by short-lived (~ 10^ yr) 
QSO activity and merging in hierarchically merging galax- 
ies are consistent with a wide range of observations in the 
redshift range 0< z <5 jHaehnelj|2o'ol . However it is not 
surprising to find a discrepancy with our analysis since we 
do not consider merging and the BH growth history that 
we find is anti-hierarchical. Our estimates agree better with 



mod els in which BH's are mainly grown by gas accretio n 
(e.g. iHaiman. Ciotti. fc Ostrikeil2003l : ICranato et al.l2004f) . 

With an estimate of AGN lifetime bas ed on BH demo- 
graph ics similar to the one presented here IYu fc Tremaind 
find r = 0.3 — 1.3 x 10* yr for the luminous quasars 



using the AGN luminosity function bv iBovle et al.l J200a. , 
The main reason our estimate is larger is that the Tued^^^l] 
(2003) luminosity function has a larger number density 
of objects, which results in longer lifetimes. Indeed, the 
quasar lifetimes co mputed using the luminosity function by 
iBovle et al.] ll2000l) are ~ 8 x lO '^ yr (dashed line in Fig.lTHb). 
in agreement with lYu fc Trem ainc (20q3)- 

Comparing to observational-based estimates, the life- 
times from this paper are in a greement with the results 
from the length of ra dio jets (see]M artinil l2003l for more de- 
tails). More recentlv. iMiller et alJ ^2003^ found from SDSS 
data that a very high fraction of galaxies host an AGN 
(~ 20 — 40%) suggesting lifetimes longer than previously 
thought (i.e. > 10* yr). Finally, our estimate is in agreement 
with the upper limit of 10® yr set by the timescale ove r which 
the quasar luminosity density rises and falls (see, e.g.. lOsmerl 
i2003). 

In summary, we estimate that local high mass BH's 
(M > 10^ Mq) have been active, on average, ~ 1.5 x 10* yr. 
On the contrary, the assembly of lower mass BH's has re- 
quired active phases lasting at least three times that much 
(~ 4.5 X 10* yr). These average lifetimes can be as large as 
~ 10® if one considers the smaller efficiency and fraction of 
Eddington luminosity which are still compatible with local 
BH's (e = 0.04, A = 0.1). 



9 CONCLUSIONS 

We have quantified the importance of mass accretion dur- 
ing AGN phases in the growth of supermassive black holes 
(BH) by comparing the mass function of black holes in the 
local universe with that expected from AGN relics, which 
are black holes grown entirely during AGN phases. 

The local BH mass function (BHMF) has been esti- 
mated by applying the well-known correlations between BH 
mass, bulge luminosity and stellar velocity dispersion to 
galaxy luminosity and velocity functions. We have found 
that different BH-galaxy correlations provide the same 
BHMF only if they hav e the same intrins i c disp ersion, con- 
firming the findings of iMarconi fc Hunti ll2003h . The den- 
sity of supermassive black holes in the local universe is 
PBH = 4.6ti ® X 10^ Mq Mpc-^ 

The relic BHMF is derived from the continuity equa- 
tion with the only assumption that AGN activity is due to 
accretion onto massive BH's and that merging is not impor- 
tant. We find that the relic BHMF at z = is generated 
mainly at z < 3 where the major part of BH's growth takes 
place. The relic BHMF at z; = is very little dependent 
on its value at Zs — 3 since the main growth of BH's took 
place at 2: < 3. Moreover, the BH growth is anti-hierarchical 
in the sense that smaller BH's {Mbh < 10*Mq) grow at 
lower redshifts [z < 1) with respect to more massive one's 
(2 ~ 1 — 3). If the correlations between BH mass and host- 
galaxy-properties hold at higher redshifts this would repre- 
sent a potential problem for hierarchical models of galaxy 
formation. 
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Figure 10. (a) BH duty cycle S{M, z) at given redshifts computed with the lUeda et al.l J2003ll luminosity function and e = 0.1, A = 1. 
(b) Average mean lifetime of acti ve BH's (AGN's) a s a function of the relic BH mass at z = computed with e = 0.1, A = 1, the 
lUeda et all i200,'j) (solid line) and iBovle et al' ] i2n0(tl (dashed line) luminosity function. The dotted line is the corresponding Salpeter 
time. The scale of the y axis on the right is the avera ge lifetime in units of the Salpeter time. The lines with the empty squares and 
triangles are the average lifetimes computed from the lUeda et alJ l200.'j) luminosity function with e = 0.04, A = 0.1, and e = 0.15, 
A = 1.6, respectively. 



Unlike previous work, we find that the BHMF of AGN 
relics is perfectly consistent with the local BHMF indicating 
the local black holes were mainly grown during AGN activ- 
ity. This agreement is obtained while satisfying, at the same 
time, the constraints imposed from the X-ray background 
both in terms of BH mass density and fraction of obscured 
AGN's. The reasons for the solution of the discrepancy at 
high masses found by other authors are the following: 

• we have taken into account the intrinsic dispersion of 
the Mbh-ct* and A^BH-ibui correlations in the determination 
of the local BHMF; 

• we have adopted the coefficients of the Mbh- 
gj, and MeH-I^bui {K band) correlations derived by 
iMarconi fc Hun3 (|2003) after considering only 'secure' BH 
masses; 

• we have derived improved bolometric corrections which 
do not take into account reprocessed IR emission in the es- 
timate of the bolometric luminosity. 

The comparison between the local and relic BHMF's 
also suggests that the merging process at low redshifts 
{z < 3) is not important in shaping the relic BHMF, and 
allows us to estimate the average radiative efficiency (e), the 
ratio between emitted and Eddington luminosity (A) and the 
average lifetime of active BH's. 

Our analysis thus suggests the following scenario: lo- 
cal black holes grew during AGN phases in which accret- 
ing matter was converted into radiation with efficiencies 
e = 0.04 — 0.16 and emitted at a fraction A = 0.1 — 1.7 of the 
Eddington luminosity. The average total lifetime of these 
active phases ranges from ~ 4.5 x 10* yr for Mbh < lO^M© 



to ~ 1.5 X 10** yr for Mbh > l(f Mq but can become as 
large as ~ 10® yr for the lowest acceptable t and A values. 
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